Spin currents, spin populations, and dielectric function 
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In Maxwellian electrodynamics, specific properties of the responses to external fields are in- 
cluded in constitutive equations. For noncentrosymmetric semiconductors, spin conductivity can 
be expressed in terms of the contribution of electric-dipole transitions between spin-split spectrum 
branches to the dielectric function. In a dissipationless regime, a spin current driven by an external 
electric field is tantamount to a background current in an equilibrium system with a reduced sym- 
metry. The importance of transients and gradients for efficient spin-current injection is emphasized. 

PACS numbers: 72.25.-b 



The growing interest in using the electron (hole) spin 
degree of freedom in semiconductor spintronics^ and 
impressing progress in the experimental study of spin- 
charge coupling through spin orbit (SO) interaction 
stimulated active research in nonequilibrium spin popula- 
tions and spin currents in semiconductor microstructures. 
The problems with achieving efficient electrical spin in- 
jection from ferromagnetic electrodes, including the effect 
of stray magnetic fields around them, stimulated devel- 
oping the concept of all-semiconductor spintronics that 
does not include any magnetic elements and is based on 
specially engineered spin injectors using SO coupling; a 
number of related ideas has been put forward^ An adia- 
batic pump for spin currents^ and a Stern-Gerlach type 
experiment with a SO beam splitter— have been reported; 
both techniques require an external magnetic field. Inde- 
pendently, an optical technique for a controlled injection 
of spin populations and spin currents has been proposed, 6 
and injection of a pure spin current (with no net charge 
current and no net spin injected) has been reported^ 
The concept of spin-polarized currents also turned out 
highly productive for metals^ and was applied to the spin- 
transfer induced switching of nanomagnets. 

Several proposals for electrical injecting spin-currents 
into nonmagnetic materials that require neither magnetic 
fields nor ferromagnetic materials have been put forward 
recently. Governale et al& and Mal'shukov et al~L cal- 
culated spin currents driven by ac electric fields, while 
Murakami et alm^ and Sinova et alw^ proposed dissipa- 
tionless spin currents driven by a dc electric field. The 
subject is of active interest, and an appropriate under- 
standing of the nature of spin currents and their poten- 
tialities for spintronics is a challenging task. 

The notion of spin current as the transport of elec- 
tron spins in a real space sounds alien to the stan- 
dard Maxwellian electrodynamics; the "magnetization 
current" c curl M has a different nature and is a part 
of the charge current. In what follows, I establish a re- 
lation between the spin conductivity and the dielectric 
function. To keep the calculations and the results as 
easy as possible, I restrict myself by free electrons. Such 
an approach allows to pinpoint the specific details of the 
band structure and specific electronic transitions that are 
responsible for spin currents, and discuss the essence of 



these currents in the framework of the standard band 
theory. It also allows to identify the electric-field driven 
dissipationless spin currentsiii^ as the equilibrium back- 
ground spin currents^ that develop in the system when 
it is subjected to a proper pyroelectric deformation. 

Because spin currents are not conserved in systems 
with SO interaction, their definition is somewhat arbi- 
trary. I apply the generally accepted and physically ap- 
pealing definition of the SC pseudotensor 
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establish its relation to the dielectric function, and dis- 
cuss the basic conditions for the generation of transport 
spin currents. Here k and v are the electron momentum 
and the velocity operator, respectively, cr is the vector 
of Pauli matrices, i,j are Cartesian coordinates, with i 
indicating the spin component and j the transport di- 
rection, and A numerates the spectrum branches. For 
T = , the integration should be performed inside the 
Fermi surface. It is an important property of the crystals 
lacking the inversion center that SO interaction splits ev- 
ery energy band into two branches. Below, calculations 
are performed for a SO-split band of two-dimensional 
(2D) electrons with a Rashba SO interaction^ A simi- 
lar approach to 3D Luttinger holes^ in a diamond type 
semiconductor will be also discussed. 

In what follows, (i) frequency dependencies of the di- 
electric function e(uj) and spin conductivity E(w) will 
be found, (ii) a relation between them established, and 
(iii) it will be shown that S(w) is directly related to the 
contribution to the real part of e(w) coming from the re- 
gion of the k -space where the lower branch is populated 
and the upper one is empty. 

The standard Hamiltonian of a 2x2 SO problem is 



H R = h 2 k 2 /2m + a(er x k ) • z, 



(2) 



with k = (k x ,k y ) the 2D momentum, z a unit vector 
perpendicular to the confinement plane, and a the SO 
coupling constant. The eigenvalues of the Hamiltonian 
Hr are e\(k) = h~uj\(k) = h 2 k 2 /2m + Xak , where A = 
±1 correspond to the upper and lower branches of the 
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spectrum, respectively; a > . The eigenspinors are 

^ k) = 7^{- iX (k x + ik y )/k)' (3) 

and the velocity operator is 

v = TT l dH R /d k =hk/m + a(ix cr)/h. (4) 

The electric current driven by an homogeneous elec- 
tric field E e~ luJt (with a vector-potential A (t) — 
e E (t) jiio ) , according to Kubo formula, is expressed 
through a retired commutator of the currents ji = evi 
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n being the 2D election concentration, and the angle 
bracket (...) indicates integration over the Fermi distri- 
bution in the (A, k ) space. Because of the isotropy of 
the problem, this integral is diagonal in the indeces . 
All matrix elements of j = ev diagonal in A cancel, 
and electric conductivity equals^ 



r(i; 



if 



H = 2^ 



<Pk 
i 



\v(k)\- 



\v(k)\-) 



1 



uj \-(k) + uj + iS lu^ (k) + lu + iS 
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with uj\\i (k) = uj\(k) — lu\i ( k ) and k± indicating 
the integration limits. When the Fermi energy is positive, 
(i > , k + and fe_ are the Fermi radii for the upper- and 
lower-branch electrons, respectively; below, all equations 
are presented for this case. Hence, this contribution to 
the conductivity comes from the interbranch transitions 
in the k -space area where the lower branch is populated 
and the upper one is empty. Using the relation 

(+| v (fe)|-) = -(-| v (fe)|+) = {ia/Tik){i X k ) (7) 

and the equation e(u>) = 4mia(uj)/u) , one comes to the 
following equations for the real and imaginary parts, 
e'(u>) and e"{ui) , of the 2D dielectric function 
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e"(uj) = ne 2 /4huj, when 2ak + < Hlu < 2ak_. (9) 

The first term in e'(cj) came from a^'(cj) and the sec- 
ond from the j -current of Eq. JSJ, to 2 — 4ire 2 n/m , 
and e' (u>) =/= only in the frequency range of the in- 
terbranch transitions. Eq. @ can be also found by the 
Golden Rule. The total oscillator strength of the inter- 
branch transition equals fso = (ma/ft 2 ) 2 /2ir . 



It is instructive to compare Eqs. JHJ and I© with the 
Kramers-Kronig relation 
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a background dielectric function that is not related to 
the electron band in question is omitted in Eq. i|l(J[l . If 
to define the Kramers-Kronig transform of Eq. @ 
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and employ the relation fc_ — fc + = 2ma/h , then Eq. 
can be rewritten as 



ul [1 - (ma/H 2 ) 2 /2 



(12) 



the ratio ma/H 2 being the characteristic SO momentum. 
To conform Eqs. (JSJ and (|1U|) one needs to supplement 
e"(ui) of Eq. © with a singular contribution 



,{uj) = {-kuj 2 p /2uj)5{uj). 



(13) 



This term reflects the oscillator strength that is hidden 
for free electrons and manifests itself in the cyclotron 
resonance, Drude absorption, etc. Eq. (12) includes the 
renormalization of this oscillator strength by the SO cou- 
pling, u 2 < uj 2 . It is reduced because a part of it was 
borrowed for the interbranch absorption of Eq. @ . 

To summarize, the function e'so(w) describes the to- 
tal contribution of the SO coupling into dielectric polar- 
izability. It includes both the direct contribution from 
the interbranch transitions near the edge of the Fermi 
distribution and the reduction of the S(oj) -part of e"(w) 
because of the oscillator strength conservation. 

Spin current J zx (t) driven by an electric field E (t) | 
y can be calculated similarly to Eq. by using J zx — 
(Hk x /m)a z . The spin conductivity equals 
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Matrix elements (A|%(fc)|A') , A ^ A', are odd with 
respect to the transposition of A and A' , while ma- 
trix elements (+|c z |— ) = (~l cr z|+) = 1 entering into 
{^\Jzx{ k )|A') are even. As a result, in a factor similar 
to the bracket of Eq. © both fractions appear with the 
same sign. Finally, the real part of T, zx (lu) equals 
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(15) 



in agreement with Schliemann and Loss^ In the low- 
frequency limit, 



E^> = 0) = e/Anh 



(16) 
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in agreement with Sinova et alH (after the difference by 
the factor h/2 in the definition of is allowed for). 

The comparison of Eq. © and (|15|l shows that the 
integrals coincide and 



8ne ma 2 /h 



(17) 



Hence, the spectrum-specific frequency dependence can- 
cels from the ratio of T,' zx {uj) and e' so (ui) . There- 
fore, spin conductivity Y*' zx (ui) , when properly normal- 
ized, acquires the meaning of the electron polarizabil- 
ity related to the transitions between the SO-split spec- 
trum branches. Remarkably, the singular low- u> part of 
e' so (ui) corresponding to £^(0; = 0) comes from the 
SO correction to w 2 , Eq. iTT^ . 

We conclude that spin currents and the SO part of the 
dipole moment, P so(t) = e so( a; ) ^ W > represent two 
aspects of the same phenomenon^ Electrically driven 
spin currents are described macro scopically through the 
SO contribution to the dielectric function. 

The implications of this observation reveal themselves 
if one takes into account that in deriving Eqs. JBJ and 
(|15|l only the interbranch matrix elements of the pertur- 
bation — (e/c)v y A y (t) were involved, while the effect of 
the electric field on the intrabranch motion dropped out. 
This is an artifact of applying Kubo formalism to free 
electrons, and doing so is an equivalent of the perturba- 
tion theory in the operator V = —eEy — ieEd/dk y in 
the 2x2 space of the spinors ip\(k) . The linear in V 
corrections to the spinors ip\(k) are 



-\(eEk x /4ak 3 )iP- X (k 



(18) 



One can now find the mean value of J zx over the new 
vacuum spanned by the spinors ip\( k ) + ip^ ( k ) . The 
contributions from the upper and lower branches cancel, 
as before, in the region k < k+ . The contribution from 
the region k + < k < fc_ equals J zx — eE/infr, in 
agreement with Eq. (|16fl . This result corroborates that 
the Kubo approach, when applied to free electrons, allows 
for the effect of the field E (t) only through the rotation 
of spinors k ) . 

Such physics corresponds to a pyroelectric deforma- 
tion along the j/-axis rather to the effect of a trans- 
port electric field E || y . This pyroelectric field lowers 
the symmetry from the group C oo V of the Hamilto- 
nian Hr to the group C s whose only nontrivial ele- 
ment is a reflection in the yz plane. There are two new 
SO invariants, cr x k y + a y k x and a z k x , in this group. 
The first can be disregarded, while the perturbation 
H' so — a z a z k x with a z rj —eE/2k 2 F results in spinors 
of Eq. H18|) and the same spin current J zx — eE/Anh 
(when fc_ — fc + <C hp ). Because a system with the 
Hamiltonian Hr + H' so is in equilibrium, these spin cur- 
rents are background (nontransport) currents. 

It follows from general arguments that time-inversion 
symmetry forbids spin accumulation in a dissipationless 



dc regime. Indeed, a linear relation between the mag- 
netization M (or spin cr ) and E is equivalent to 
a magneto-electric effect that is generally forbidden 
because E is real while M is imaginary with re- 
spect to t -inversion. For the magneto-electric effect 
to exist, this symmetry should be violated by a proper 
magnetic structure^ electron scattering, 21 or a finite 
frequency, w— Remarkably, spin polarization of free elec- 
trons that develops in the direction (E x z) || x and 
diverges as w" 1 when u> — > , is cut-off by the momen- 
tum relaxation time t v ^ this effect has been observed 
recently^ 

In noncentrosymmctric crystals, the existence of equi- 
librium background spin currentgi&i^ is compatible with 
the t — > —t symmetry. In fact, this symmetry re- 
quires that momenta and spins be reversed simultane- 
ously, hence, it does not require that the currents of 
the particles with a given spin vanish^* Therefore, in 
a dissipationless regime Kubo formalism maps the real 
system driven by an electric field E onto an auxiliary 
equilibrium system of a lower symmetry. Both systems 
are described by identical equations because intrabranch 
dynamics has been eliminated. Spin currents flowing in 
the mimicking system are background currents. The op- 
tical analogy suggests that they set an upper bond for 
dc transport spin currents in the real system because the 
substitution u> — > — iT in Eq. JTSJ, T being a proper 
decay constant, results in a decrease of T,' ZX (T) with T . 

The mapping of the real system onto an auxiliary equi- 
librium system is helpful because the problem of dissi- 
pationless currents cannot be posed rigorously in the ab- 
sence of a strong magnetic field, and the mapping clarifies 
the assumptions underlying it. 

Background spin currents do exist because the oper- 
ator Jij is real with respect to t -inversion. For the 
same reason E^ -(w) is related to e'(u>) , the dispersive 
part of e(uj) , rather than to its dissipative part e"(w) . 
These currents are a reality. Nevertheless, they "do not 
work" as spin sources because a pyroelectric element of 
a circuit cannot inject spins at equilibrium. To realize 
how the currents of this sort can be put to work, it is 
instructive to consider a simple classical analogy. The 
momentum flow 11^ = HiVi{l)pi{l) is also real with re- 
spect to t -inversion, I numerates particles inside a unit 
volume. For an equilibrium gas IIm = P , the pressure. 
In macroscopic terms, the equation for the momentum 
flow of an ideal fluid is d(pvi)/dt = —dUij/dxj , where 
n,y = PSij + pviVj is the tensor of the momentum flow 
density, p is a density^ It is a gradient of Ily that 
produces an acceleration, a flow, and transients. A simi- 
lar approach is valid for spin currents. E.g., Mal'shukov 
et aim have shown that modulating a = a{t) results 
in injecting spin currents. In a diffusive regime they are 
controlled by the ratio ui/T i& For abrupt changes in a , 
ballistic pulses can be anticipated. In Stevens et ali ex- 
periments, spin-current pulses spread from a small spot 
where they were generated by a laser pump. 

Dynamics and propagation of spin populations and 
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spin currents should be based on a theory including dis- 
sipation, and different approaches to this problem have 
already been advanced»ifl»i2*2£i2L22i In particular, such a 
theory should provide generalizations of Eq. I|17|) . When 
SO coupling is strong, akp ~ /i , the only characteristic 
time is r p , the momentum relaxation time. Long relax- 
ation times ts of spin populations^ can be achieved 
when SO coupling is weak, akp <C fJ. , and t p is short, 
t p <C Ti/ akp . Then T5 , controlled by the Dyakonov- 
Perel process^ is long, 1 w T p (2akp /K) 2 . In this 
regime, a quasiequilibrium in orbital degrees of freedom is 
established, 30,31 and spin dynamics in an external electric 
field is strongly influenced by momentum scattering^ In 
the absence of spin populations, spin currents decay at 
the t p scale. Some of the emerging problems are similar 
to those of the physics of spin photocurrents^ 

Holes in centrosymmetric crystals of the diamond type, 
described by a Luttinger Hamiltoniani^ 

Hl = [(71 + 5 72 /2)fc 2 - 2 72 ( J ■ k ) 2 ]/2m, (19) 

are in some aspects similar to the electrons of the Hamil- 
tonian Hr . Their spectrum consists of two (twice de- 
generate) branches known as heavy and light holes with 



the energies (71 =p ^l2)k 2 /2m . The momentum J is 
described by 4 x 4 matrices of the angular momentum 
J = 3/2 . The currents Jij defined similarly to Eq. (JJJ 
with <r > J have the meaning of angular momentum 
currents and are related to interbranch transitions, in 
this case from the heavy hole to the light hole branch^ 
Therefore, most of the above conclusions are applicable 
to this system. The main difference is related to the 
fact that for Luttinger holes the notion of "spin" should 
be understood generally, as a total angular momentum 
rather than the physical spin, the constant 72 is only 
weakly influenced by the physical SO coupling, 16 and the 
heavy-light hole relaxation time is never long. 

In conclusion, spin currents in noncentrosymmetric 
semiconductors are accompanied by electric dipoles re- 
lated to electronic transitions between spin-split spec- 
trum branches, and the polarization of eigenspinors by 
an electric field produces background spin currents. From 
this standpoint, the importance of transients and gradi- 
ents for efficient electrical spin-current injection has been 
clarified. 
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